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Abstract 

We have estimated higher order quantum gravity corrections to de Sit- 
ter spacetime. Our results suggest that, while the classical spacetime metric 
may be distorted by the graviton self-interactions, the corrections are rela- 
tively weaker than previously thought, possibly growing like a power rather 
than exponentially in time. 
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As is very well known Einstein field equations permit an addition of an extra term 
proportional to the metric tensor, g^ u : 

R^u - 7,9tMvR = SttGT^ + Ag^ (1) 

Here T^ u is the energy-momentum tensor of matter and A is the so called cosmological 
constant. (We put an bar over all quantities here to distinguish them from the rescaled 
ones, see Eqs. (fUD). The last term in Eq. ([!]) may be interpreted as the vacuum energy- 
momentum tensor. Any natural estimate gives A or p vac = A/8irG much larger (by 50-100 
orders of magnitude) than the observed upper bound in the present day Universe. This 
mysterious discrepancy is known as the cosmological constant problem and presents one of 
the most interesting challenges in the modern physics (for a review see Refs. [TO, 0]). 

de Sitter spacetime is a solution of the Einstein equations with a dominant cosmological 
term. Written in the Robertson- Walker form for the special case of the spatially flat section 
the metric takes the form 

ds 2 = dt 2 - a 2 (t)dr 2 (2) 

where a{t) = exp(Ht) with H = WA/3. The assumption of spatial flatness is not essential 
and the results obtained below are true also for open and closed geometries. It is convenient 
to rewrite the metric in terms of conformal coordinates where, up to an overall scale factor, 
it has the Minkowski form: 

ds 2 = a 2 (r)(dT 2 -df 2 ) (3) 

The conformal time r is related to the physical one as dr = — exp(—Ht)dt and a(r) = 
— 1/Ht. Note that when t tends to future infinity, r tends to 0—. The de Sitter metric 
possesses the same degree of high symmetry as Minkowski one and for this reason, it is the 
simplest non-trivial curved background for a quantum field theory. Despite its simplicity, it 
has several very interesting properties and in particular generates an infrared instability of 
a massless scalar field (f> minimally coupled to gravity. It was shown || || that the vacuum 
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expectation value (VEV) (0 2 ) is singular at the zero mass limit, (0 2 ) ~ H 4 /m 2 . With 
the advent of the inflationary scenario, this phenomenon was rediscovered in a number of 
papers |5|, || [jj where it was argued that (0 2 ) rises linearly with time in the zero mass 
limit. This quantity however does not have a direct physical meaning and, in particular, the 
contribution of the field <p to the energy density does not rise with time. Nevertheless, this 
infrared instability afflicts the scalar field propagator and needs to be resolved to make a 
sensible quantum theory. Note in passing that this instability is a result of broken conformal 
invariance. To ensure the latter one needs the coupling to the curvature scalar, R(p 2 /6, 
which gives rise to an infrared cut-off. In fact, any nonzero mass or coupling £i?0 2 with a 
nonnegative coefficient is sufficient for infrared stability. In any case, this infrared divergence 
is a sign that the true vacuum state is not de Sitter invariant. 

Conformal invariance is also broken for gravitons || and, for this reason, there should 
be infrared instability of the de Sitter vacuum due to quantum gravity effects. If this is 
indeed the case, the solution of the long standing cosmological constant problem may be 
found in this direction. One-loop graviton quantum corrections to the de Sitter metric were 
considered in this connection in Ref. || where it was found that they are time independent. 
However, it was suggested that higher loops may show evidence of this infrared instability. 
Indeed, it was claimed recently in several papers |XD| , |TT| , |T2| ] that higher loop effects are 
much stronger, giving corrections that rise exponentially with physical time t or as a power 
of 1/r with conformal time. This is a very exciting result and, if confirmed, would mean 
that de Sitter space cannot exist indefinitely, opening a beautiful way for the solution of the 
cosmological constant problem in the framework of the normal quantum gravity without any 
drastic assumptions. 

Here we have reconsidered results of papers [0, [IT], [12] using a different formalism and 
have found, unfortunately, that such a strong instability does not set in. The corrections at 
most behave as a powerf] of log |r| ~ t. 

2 We understand that these authors now also do not get power law singularity but only powers of log r. (R. 
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For the conformally flat background metric we introduce the quantum graviton field 
in the usual way 

9tw = a 2 (i~)gv,v = a 2 (r)(?7^ + h^). (4) 

The issue is whether, taking into account loop corrections, the VEV (h^) is non-zero and, in 
particular, whether it is time-dependent. This would suggest that the background de Sitter 
metric is not self-consistent, although it must then be shown that this is a physical effect by 
calculating, for example, the curvature for the modified metric. Even then, it is important to 
ascertain whether the inconsistency involves strong coupling or whether it is an essentially 
negligible effect. The Einstein action (with nonzero cosmological constant) in terms of the 
new metric g^ can be rewritten as 

11/3 

A — — I d 4 x^(R + 2A) = — / d 4 xa 2 ^(R + 6^1 + 2Aa 2 ), (5) 
with k 2 = 16n/m 2 Pl . This implies the following equations of motion: 



R P v - ^9^R - Aa 9/au H ~ 2 — + 9/u, — ! — ) = (6) 



where tensor indices are raised and lowered with the new metric g^ v . To zeroth order in 

h^, we get the usual equation for the scale factor of the classical background metric, with 

solution a(r) = —1/Hr (see notation after Eq. (Q).) 

It is of course well-known that this is a non-renormalizable quantum field theory, but 

we regard the Einstein action as the first two terms in an infinite series of local operators 

of increasing dimension. This is an effective field theory, presumed valid on scales below 

mpi, and is renormalizable in the sense that all divergences involving these vertices may 

be absorbed in a renormalization of one of the infinite number of coefficients of these local 

operators. While we have not written such terms explicitly, they should be understood to be 

present. In the following, we will deal exclusively with renormalized couplings and operators, 

with the tacit understanding that the counterterms are included in our interaction. To 
Woodard, private communication) 



consistently quantize the theory, we must add a gauge-fixing term C g f and the corresponding 
Faddeev-Popov ghosts Cfp- We choose C g f = — \F ii F v Tf u as in Refs. [H], [Tl|, with 



1, _n,„a, 
-/ 
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F, = a{r){h^--h, + 28lh^) (7) 



The ghost Lagrangian Cpp may be found in Ref. [pL2|| . Here and subsequently when we 
consider perturbation theory in h^, the indices are raised with the Minkowski tensor 77^. 
With these gauge conditions, the linear part of the equation of motion for h^ u has very simple 
form: 

2 2 2 

Ku,a V,0 + -o ^ + 5° V) oV^OO = ( 8 ) 



One can easily verify that time components ho^ (/i = 0, 1, 2, 3) are conformally invariant in 
the sense that the reseated functions xo^ — a ( T )hon satisfy the free field equations of motion, 
d 2 Xo, = 0. 

More interesting are the space-space components of metric hij. The equations of motion 
are not diagonal for them but after a simple linear redefinition = — 5ijhoo they are 
diagonalized and have the form: 

fia ~ -fiifi = (9) 
T 

This is the same equation which is satisfied by a massless minimally coupled scalar field. 
The solutions with definite momenta can be written as 

0fe(r, x) — H — — ^ exp(ikx — ikr) (10) 

These nonconformal modes differ from the conformal modes by the presence of the i/k term. 

The conventionally normalized Heisenberg field operator for quantum fluctuations of the 
gravitational field is defined as ip^ = h^ v j ' k. This may be expressed in the standard way in 
terms of creation-annihilation operators a^k) by the decomposition:^ 



J See, e.g., Ref. fllq] for further details. 



with a and a' satisfying the commutation relations 

a /J , L ,(k),a' flv (q) = (2ir) 3 8 3 (k — q) (not summed on /i, v) (12) 

All other commutators are zero. For the nonconformal modes, the wave function Tpk(r,x) is 
given by Eq. (|10D , while for the conformal modes, they are plane waves Hr exp(ikx — ikr). 
Respectively, the propagator for a conformal field (for example for the field ho a , see above) 
is proportional to that in flat space. It means that the evolution of ho a is trivial and 
in particular the corresponding Green's functions are obtained from the ones for the flat 
spacetime by the rescaling: 

G dS ( x , x') = , } ' J = H 2 rr'G flat (x, x'). (13) 
a{T)a(T') 

Now, propagators (two-point functions) for nonconformal fields exhibit new features since 
the solutions ([TOP blow up at small k as compared to the solution for conformal fields. In 
particular, this implies that the anticommutator correlator G\ is in fact not defined in the 
infrared. To be explicit, it is given by 

Gi(r, r', k) = ^-[{tt' + i cos(A;Ar) - ^ sin(fcAr)], (14) 



where Ar = r — r'. If we calculate the Fourier transform of fll4| ) then it is logarithmically 
divergent at small k and undefined. However, derivatives of it are well defined. 

It might also worth mentioning that in a de Sitter background, the Green function for 
nonconformal field can be also given in a closed form in four dimensional notations. Namely 
we have 

A , f d A k exp i(k,x — x') f d 4 k exp iik.x — x') . . 

= TT J — p J — ^ — (15) 

where, as usual, depending on the ie prescription we get either Feynman or retarded, or 
advanced propagator. (The derivation of ( Jl5| ) follows the lines presented in the appendix to 
Ref |L7| where the first terms in the WKB expansion are considered for arbitrary background 
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field. In a de Sitter background, a similar argument fixes the exact propagator.) In particular 



the Feynman propagator takes the form[|Tl 

Gf(x, x') 



An 2 



TT ' 1 I / A2 

m(x — x ) 



(16) 



(x - x') 2 2 

Note that the argument of the log is not defined, as a manifestation of the same infrared 
instability. What is specific about the de Sitter background (and £ = 1/6 or 0) is that the 
expansion in k~ 2 in Eq. ( |15D terminates on the second term. 

In fact, this infrared blowup of the Feynman propagator is the basic observation which 
gave rise to the hopes that quantum corrections to the de Sitter metric grow with time. 
Indeed, the quantum propagators are not vanishing in causally forbidden regions and may 
bring information on the overall expansion in the future, characteristic of the de Sitter 
solution, revealing in this way an instability of pure classical solution. We are going to 
scrutinize this suggestion. 

The best way to approach the problem is to keep as close to classical consideration 
as possible. Indeed, the classical solutions are known to be stable against perturbations 
| |18| |. The crucial difference between classical and quantum problems is that development 
of classical fluctuations is governed by the retarded Green function Gr which is free of the 
infrared divergencies mentioned above. The Green function is obtained in the standard way: 

Gr(x,x') =i6(T-T')MT,x)MT',2')]) (17) 

where we have suppressed the various tensor indices. Although commonly expressed in terms 
of a VEV, it is important to recognize that Gr for the linearized theory is a purely classical 
construct that may be obtained directly from the classical equations of motion. Accordingly, 
it is completely independent of the definition of the vacuum state. 



A simple calculation gives [12 



G R {x,x') = H 2 [rr'G R lat (x - x') + #(At)0(At - r)/4n] 
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where r = \x — x'\, At = (r — r') and 

G% at = 6(AT)5(r - At) /4vrr (19) 

is the retarded Green's function in the flat spacetime. We suppressed here the evident tensor 
indices. The second term in Eq. fll8f) is connected with the broken conformal invariance. For 
what follows it is essential that it vanishes for At = 0. In the mixed (r, k) representation 
we find 

G r (t, t', k) = H2() ^ T) 1(tt' + 1) sin^Ar) - ^ cos(*Ar)], (20) 

and in the limit if A; — ^ the retarded Green function for a nonconformal field is no more 
singular than for a conformal (or in flat space). 

The Heisenberg picture of quantum field theory is most close to the classical one and 
we expect that the infrared problem is most easily treated within this approach. With the 
initial condition Vv^o?^) = Vu™ ( r 0; the operator equation (||) can be rewritten in the 
integral form: 

^V(r, x) = ^(r, x) + — f dr' I d'x'Gf^r, x; t', ^(H^V^, x') (21) 

mpi Jt j 

where V a p{T, x) is the interaction term. More specifically, 

v, v = ^fvw + - iM%, + 2T- U r% - 2 r^ri + 2 V - 

(22) 

where r°„ = + ~ i J, / Z)/^ * s the Christoffel symbol associated with the metric 

9fiu = Vfiv + K Vv ^° ^ rs ^ or der in Note that the last two terms that are explicitly 

singular in r are multiplied by combinations of fields whose wave functions, in lowest order, 
vanish when r — > 0. This result persists in higher order, as we shall argue below. 

We will be looking into the perturbative solution of ([H]). The use of fl2*l|) insures that 



at least one propagator attached to each vertex is a retarded one and is smooth in infrared 



We remind the reader that counterterms are implicitly included in V^. 



limit. Another important advantage of the Heisenberg representation is that the state vectors 
in Hilbert space remain constant and the evolution of the system is contained in the time 
dependence of the field operators. Therefore one may study the time evolution of the physical 
quantities such as, e.g., energy density as the expectation values of the corresponding field 
operators over the initial quantum state \in) which we take to be a no-particle, de Sitter 
invariant state. The time evolution of this state is quite a complicated problem by itself 
since the notion of particle production is ambiguous in general relatively (for discussion see, 
e.g., flT5|). In the Heisenberg picture we do not confront this problem. 

Although analytically quite complicated, it is straightforward in principle to develop the 
perturbative solution of the integral equation ([H]) (See Fig. la.) Each three-point vertex 
carries a factor of 1/mpi, so, for consistency, we should also include in higher-point 
vertices obtained from the expansion of Eq. (|6|) as well as vertices of the same order from 
higher dimensional operators implicitly included in the action Eq. (^). We will discuss 
that in more detail in a subsequent paper where we will also consider in greater detail 
the renormalization of the theory. Here we concentrate only on a possible existence of 
singularities at r = 0. The contribution to ip from the conformal Green's function reads 

t r dr' 

A con / = / — - 

mpi J to r 

where \r\ = t — t' and the angular integration is made over directions of f. The nonconformal 
part coming from ( |18|) contains two-terms, the first of the form of Eq. (|2~3"D and the second 
given by 

I f T (It 1 r( T ~ T ') 
^nonconf = -. / drr 2 'dilV '(/, X + r) (24) 



£^-(r-r') I d{lV(r',x + r) (23) 



Our interest is in the VEV of this expression. Using translation invariance, one may simply 
set the spatial arguments to to obtain 

l pi J TO 



1 r T dT 1 

(A„/) = / — (r-r') 3 (VV,0)> (25) 

3m pi J to t 1 
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The generic form of the loop expansion is illustrated in Fig. lb. In calculating the VEV, one 
encounters quantum correlation functions of the type 

Gi(t,t',x - x') = (Q,in\{i> in (T,fy,i> in (T',x*)}\0,in). (26) 

These correlators are not time-ordered and replace the Feynman propagators of the more 
familiar in-out formalism. Unlike the retarded propagator or commutator function, these 
are truly quantum-mechanical amplitudes that do not vanish for spacelike separated points. 
They depend upon the definition of the vacuum state |0, in), about which we shall have more 
to say shortly. 

The dimension of is mass 4 . At one loop, it involves a single graviton propagator, 
which is proportional to H 2 . It can be shown that each additional loop brings out another 
factor of H 2 /m 2 Pl . Thus there are two powers of mass to be accounted for. Since the only 
remaining scale on which the VEV (0, m|V^ v (r, 0)|0, in) can depend is r, one would therefore 
naively expect it to behave as 

H 2 {H 2 \ L - 1 , s 

(27) 



m 



pi 



times possible powers of log|r|. This behavior, when inserted into the integral equations, 
implies that the conformal modes are finite as r — > 0, while the nonconformal modes may 
diverge as a power of log \r\. This would suggest that the de Sitter background would be un- 
stable to quantum fluctuations, since the fluctuations would be growing with time. Although 



this is a much weaker singularity at r = than previously suggested, [11], |I2"| , this may 
nevertheless lead to a breakdown a perturbation theory and leave open the question of the 
ultimate future of the de Sitter metric. 

It might worth emphasizing that it is not only the use of Gr that softens the infrared 
behaviour but the form of as well. With one exception, the nonconformal fields ipij enter 
the expression V^ v in the form ipij t(X , and the derivative renders their contributions infrared 
finite. The one exception is the second-to- last term in Eq. (|22|), which includes ^T^ /t. 
The connection T?- involves space derivatives of conformal modes ipoi as well as the time 
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derivative of the nonconformal modes ipijfi- Despite what one might think, if>ij t o can be 
shown to satisfy a conformal field equation, so that all these terms in T^- involve an explicit 
factor of r and satisfy the conformal integral equation. Thus, this term in is a kind of 
cross product of a nonconformal field with a conformal field. It is only the expectation value 
of squares of conformal fields that manifest infrared problems, so one would not expect this 
to give trouble. In conclusion, there are no infrared divergences in quantities of interest to us 
here, and we may ignore for our purposes the fact that the true vacuum state is not de Sitter 
invariant. 

Let us consider higher loop corrections to (V(r, 0)). The two-loop corrections take the 
generic form 

1 r' dr\ r' dr 2 f d 3 k d 3 q 
m^Jro If Jr If J (2tt) 3 (2tt) 3 

x d 2 Gx(ri, r 2 , qj&dfa, r 2 ,k - q) . (28) 



8) = i £' f I' f / W, n, S) r 2 , -k) 



Here we have gone over to a mixed time-momentum representation, which is more convenient, 

since momentum is conserved, and power-counting is more familiar in momentum space. Gr 

is, as before, the retarded Green function, while G\ is the anticommutator functionP] given in 

Eq. (|26l) . Each partial derivative in Eq. (|28|) is to be interpreted either as a derivative with 

respect to time or multiplication by one power of momentum. As mentioned previously, these 

remove the potential infrared problem normally associated with the presence of G\. With 

regard to the momentum integration, there will be various terms. In a mass-independent 

renormalization prescription, such as dimensional regularization with minimal subtraction, 

all divergences proportional to a power of the cutoff are completely removed by the countert- 

erms. Thus, we need only concern ourselves with terms that are logarithmically divergent. 

Such contributions cannot change the power behavior of the time dependence, so that we 

simply have the result that one would guess on the basis of dimensional analysis. Thus, at 

worst, F 2 ~ (log \t'\) p /t' 2 for some power p. When inserted into the integral equation for 
There are also two-loop contributions from the commutator function. 
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the nonconformal iffy, this yields a result that behaves at worst as (log |t|) p+1 . 

Thus it is safe to say that quantum correction may bring only logarithmic dependence 
on t. Terms proportional to powers of In |t| are not inherently problems for perturbation 
theory. Consider, for example, the effect of a finite renormalization of the curvature, so that 
the scale factor a(t) might be of the form exp[(H + 5H)t]. Expanding in powers of of the 
perturbation, gives 

exp[(H + 5H)t] = exp(Ht)[l + SHt + i^L + ....] (29) 

In terms of the original conformal time Ht = — ln(iJ|r|), these corrections take the form 
of a power of ln|r|. Thus, it may be that logarithms arising from loop corrections, while 
suggesting a breakdown of perturbation theory, are merely finite renormalizations of the 
Hubble constant. To understand whether logarithmic corrections are a true instability rather 
than simply a finite renormalization, one would have to show that the breakdown is not 
simply due to large logs that can be summed up (as is often done using the renormalization 
group). Their physical interpretation may be simplified by computing their contribution to 
a gauge invariant quantity, such as the curvature. 

Let us now outline an argument according to which these logs may in fact severely 
constrained further. To this end, we return to the consideration of two-loop correction ( p8| ) 
and Fig. 2. We have not completely settled the issue of how many logs are present, but we 
can present an argument that no logs arise from the quantum loop integral. To this end, 
we will show that the one-loop corrections to the effective action does not contain a In |t| 
term in the two-point function that enters the two- loop calculation. The r.h.s. of Eq. (p8|) 
indicates in fact two distinct steps in the calculation. First, we calculate the one- loop effective 
action which involves acausal propagators. Then we iterate this action with the help of the 
retarded (i.e., classical) propagators. Our procedure so far was to look for log dependence 
of the generic term fl2"8p. We did not check that the result adds up to a general covariant 
function. In fact, as we will demonstrate now, this constraint may eliminate some infrared 
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logs. 

To this end let us consider the effective action generated by a scalar field of mass m and 
parameter £ (see, e.g., chapter 6 of book ]TB|): 

L di v = - (4tt)-7 (—- - + i( 7 + ln(m7// 2 )) • a 2 (x) (30) 

where 7 is Euler's constant, and 

a 2 (x) = ^R a ^sR a ^ 5 - ^R a pR aP -l^- m aa + l(l- 2 R 2 - (31) 

Here R a p 1 s ) R a /3, R are the Riemann tensor, Ricci tensor and scalar curvature, respectively. 
Since quantum part of the gravitational field, h^, is a collection of fields with £ = or £ = 
1/6 (see above for a discussion,) Eq. (^TJ) can be directly used to evaluate the effective action 
generated by quantum gravity. Note also that to regularize the expression in ultraviolet and 
in infrared one introduces both dimensional regularization, n^4, and a nonvanishing mass, 
m^fl. 

Eq (|30|) is valid in arbitrary background field. For our purposes we need to calculate 
Ra/3-yS for the metric specified by (|]) and expand the result in ip^. What is specific about 
de Sitter background is that 

This means that in the limit n — > 4 only anomalous terms survive where, by the anomalous 
terms, we now understand terms proportional to [n — 4) _1 . Thus the correct result for the 
variation is not zero but local terms. 

This kind of technique is widely used to fix the trace anomaly (see, e.g., ||16|| ) in arbitrary 
background. What is specific about the de Sitter background is that there are no "normal" 
terms in one-loop expectation value of the energy-momentum tensor (T^); the whole of 
{Tpv) is anomalous both for £ = 1/6 and for £ = 0. While considering logs piecewise, 
without collecting terms into R 2 or other invariants, we lose the property (|32|) which in fact 
extends to second derivative as well, as we argue next. 
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Indeed, what is special about the terms in the r.h.s. of eq ([31]) is that, upon multiplication 
by V~9 they do not depend on the Hubble constant H. Therefore 

5" 



■(VgR 2 ) = o (33) 



(shy 

for any n. The same is true for all the invariants quadratic in R a p~f& entering 02 (x). On 
the other hand, in conformal coordinates, the variation of H is equivalent to a variation of 
metric proportional to the metric of de Sitter space itself: 

5 2 
jjj9^u = -jjQuv (34) 

Thus there exists a direction along which even finite variation of leaves y/—ga,2(x) in- 
variant. High symmetry of the de Sitter solution implies then that variations of y/— ga,2(x) 
are strongly constrained for arbitrary change of g^ as well. In particular, the first variation 
is always proportional to g^ v itself: 

S(V=9°*(*)) Aimer =r(x)<rs gia , (35) 



de Sitter 

where r(x) is function of coordinates. Since we know that the variation vanishes for the 
particular choice (|3~4D of Sg^ we conclude that in fact r(x) = 0. 

For our purposes here we need the second variation of the effective action above. This 
vanishes again for n = 4. Indeed the second derivative generically takes the form: 

S 2 (y^a.ix))^ smer = N^SSg a0 g jS , (36) 

where N a p r s is a tensor which in case of de Sitter space is constructed on g^ alone. Moreover 
it vanishes upon multiplying by the variation of the special form (B3J). It means, 



N a f3~ ( 6 ~ 9-yagsp ~ 9^9 So. (37) 

and the second derivative ( |36"D vanishes for any variation Sg^ symmetric in \xv. Thus we 
have shown that the second variation of the action can produce only anomalous local terms. 
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Such terms are to be included in renormalization of the bare action and are not relevant to 
infrared logs. So far we relied on infrared regularization by finite mass. If we put the mass 
equal to zero from the very beginning, then one may argue [Tj | that one has ln(i?) instead 



of ln(m). This form of the infrared cutoff does not modify our conclusions. Indeed, if one 
calculates variations of such an effective action, the result is again local terms. 

To summarize, we have demonstrated that quantum corrections to the classical de Sitter 
solution generated by higher loops in quantum gravity can be at worst powers of logs of the 
conformal time r. Moreover calculation of these logs is no longer a pure infrared problem. 
We presented then further arguments based on consideration of the one-loop effective action 
according to which one-loop effects in de Sitter background reduce to renormalization of the 
action for ip^ V) with no logs involved. One needs also to include the effects of the FP ghosts 
as well as higher-point vertices of the same order in l/rripi. 
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his results prior to publication. We are also grateful to A.I. Vainshtein and M.B. Voloshin 
for useful discussions. A. Dolgov is grateful to the Department of Physics of the University 
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Figure Caption 

Figla: Perturbative solution of Eq. 21. The line denoted by R corresponds to the retarded 
propagator. 

Figlb: Loop corrections to metric fluctuations tp^. 
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